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Abstract: The consequences of the sudden change in the coupling constants (quenches) 
on the phase structure of the theory at late times are explored. We study in detail the 
three dimensional 6 model in the large N limit, and show that the 6 coupling enjoys a 
widened range of stability compared to the static scenario. Moreover, a new massive phase 
emerges, which for sufficiently large coupling becomes the dominant vacuum. We argue 
that these novel phenomena cannot be described by a simple thermalization effect or the 
emergence of a single effective temperature. 
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1 Introduction 

Unlocking the problem of out-of-equilibrium dynamics of a quantum coherent system is one 
of the fundamental questions in quantum phyiscs. This is particularly true in the context of 
quantum field theories, where many important questions have so far been addressed mainly 
in a static scenario, such as the renormalization group flow, the phase structure of vacuua 
and critical points. There are also interesting questions that spring directly from a non- 
equilibrium system, such as the mechanisms of relaxation, the time-scale over which this 
occurs, and the existence of an effective description at late times. These problems, given 
their fundamental role in field theories, naturally appear in many places. For example, 
it is not surprising that non-equilibrium dynamics are important in cosmology, in the 
evolution of the early universe. The RHIC experiments, involving the relaxation of the 
quark-gluon plasma, is another such example. Dynamical systems appear frequently in the 
context of condensed matter physics. Recently, the study has been rendered particularly 
pertinent experimentally due to new advances in the control of cold atomic gases [1-4]. 
For the first time, we are able to observe minute details of the evolution of a system 
that retains its quantum coherence for sufficiently long periods of time. One class of 
situations that has been subjected to intensive studies is called the quantum quench, in 
which a particular external field, or parameter, of the system is changed abruptly. An 
example is a sudden change of the external magnetic field to which the atoms couple. 
These experiments have inspired a flurry of theoretical activities, most notably initiated 
by Calabrese and Cardy [5]. Previous works however, have been concentrated on free 
field theories, one-dimensional interacting theories and integrable models, see e.g. [6-12] 
Attempts to understand interacting theories in higher dimensions by considering the large 
N limit of a </> 4 theory have been made in [13]. (See also the related problem [14]). 

Previous studies of the quenches concerned mostly about the relaxation of the system. 
One central issue is whether the system thermalizes and therefore describable by an effective 
temperature at late times. It is however an open problem if thermalization occurs at all, 
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and if it does not, which is shown to be the case in many integrable models and in some 
cases even interacting models (see for example [10]), whether there are convenient effective 
descriptions of such systems and observables or effective parameters that characterize their 
behavior. This leads us to the current investigation of the phase structure of some out-of- 
equilibrium state, which in the scenario concerned is prepared by a quantum quench. This 
should be contrasted with the usual notion of the phase structure of a given Hamiltonian, 
which is a property of its ground state. Here, we have to deal with a state, which, while 
settling to some static equilibrium in the far future, does not resemble a thermal state, nor 
is it able to relax to the ground state because of its isolation and energy conservation after 
the quench. It is therefore only natural to consider fluctuations about such a special state 
as oppose to the ground state, and determine its corresponding phase structure. 

In this work, we demonstrate that this phase structure differs significantly from that 
of the ground state even at late times as the system approaches equilibrium again. 

In particular, we explore the gQ(j) e theory at the tricritical point, i.e., when all dimen- 
sionful couplings except the physical mass immediately after the quench are tuned to zero. 
What is special about this model is that it was shown to possess an ultraviolet fixed point 
g$ = 192 using the 1/N expansion about N = oo [15, 16]. This fixed point however lies in 
the instability region of the model where the non-perturbative effects dominate 1 [18]. The 
latter implies that the theory is always driven into the unstable region by the /3-function 
and therefore apparently does not make physical sense. We revisit this theory in the con- 
text of quantum quenches. To that end, we employ the methods introduced in [13], where 
the effect of a quench is incorporated as a boundary condition on the fields. Assuming 
that the system does settle down, we then self-consistently compute the effective potential 
which defines the notion of phase structure of the theory at late times. We thus obtain a 
corresponding phase diagram, which surprisingly is modified dramatically in comparison 
to the unquenched case. The region of stability is substantially widened such that the UV 
fixed point of the /3-function now lies well within. Moreover, a new stable minimum in 
the effective potential emerges when the coupling constant exceeds the upper bound of the 
stability range in the static theory. The new vacuum starts life as a meta-stable phase, 
but then becomes dominant for sufficiently large values of the coupling. In particular, the 
effective mass of the new phase increases as the coupling increases, and eventually diverges 
when the coupling hits the boundary of a newly established range of stability. 

In the following, we will employ large iV expansion and study in detail the quench 
dynamics of the massive scalar 0(N) vector model. 

2 Quenching the scalar model 

The scalar O(N) vector model consists an iV-component scalar field <fi. For simplicity we 
assume that initially the theory is free and the system is prepared in the ground state of a 
free hamiltonian |^o)- Ait = the marginal </> 6 interaction as well as relevant cj) 4 interaction 
are instantaneously switched on, and at the same instant the bare mass parameter of the 

1 See also [17] for recent analysis of the /3-function in the case of three dimensional Chern-Simons theories 
coupled to a scalar field in the fundamental representation. 
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field jumps from po to p. The action of the system after the quench is given by 

S(<f>) = \ J d* x [d,4>d^ - M V - |^ 2 ) 2 - ^- 2 i4> 2 ?] ■ (2.1) 

Since parameters of the theory are changed abruptly rather than adiabatically, one 
needs to resort to the well-known Keldysh-Schwinger, or in-in, formalism for non-equilibrium 
quantum systems. In this formalism the integration over time coordinate t in the path in- 
tegral starts from some initial time ti, extends to some final time tf and then goes back to 
ij. Correlation functions are path ordered. In this approach one needs to impose boundary 
conditions at t = ti. In our case we impose that the initial state at ti = is given by l^o)- 

The expectation value of an arbitrary operator 0(t) is given by 

(*o|0(*)|*o> = / D4>6(t)e iS ^ , (2.2) 

JCTP 

where for brevity we use the following notation to designate the closed-time-path (CTP) 
integral measure 

/ D<j> = I DfaV (fa) [ Dfa%{fa) I** D<f>, (2.3) 

JCTP J J J<j> t 

where (pi and fa denote the values of the scalar field 4> at the end points of the time contour, 
whereas ^o(fa) = (^il^o) an d similarly for the complex conjugate ^^(fa). 
Introducing the following identity into the path integral 2 

I~ / Dp5{(p 2 - Np) ~ f DpDX e -U^x\(^-Np) ? ( 2A j 

JCTP JCTP 

(tf |0(t)|*o}= / D(j) I DpD\6{t)e iS ^ X) , (2.5) 

JCTP J 



yields 



where 



S(fa P,X) = \J d 3 x [d^<f> - (p 2 + A)0 2 - N 9 -^- - N 9 ^- + NpX] . (2.6) 
Performing now the Gaussian integral over <f> leads to 

(*o|0(t)|* > = / DpDX6{t)e iNS *ff^ x) , (2.7) 

JCTP 



with 

^P 94 2 96 3 



S eff (p,X)= / d 



+ -Trlog(D + / u 2 + A) . (2i 



2 We keep CTP label in the path integral over p and A to emphasize that the delta-function is inserted 
at each point of the Keldysh-Schwinger contour. Obviously there are no boundary conditions associated 
with p and A. 
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The first thing to note about the above expression is that boundary conditions are 
now encoded in the functional trace. Secondly, this trace explicitly depends on the integra- 
tion parameter A, and this in turn renders evaluation of the remaining path integral very 
difficult. 

However, in the limit when N is large while 54 and g& are fixed, the right hand side 
of (2.7) is dominated by the field configurations which minimize (2.8), i.e., solutions of 
the corresponding classical equations of motion. The effective mass can thus be evaluated. 
This is often called the stationary phase approximation. This gives 

m l = fJ 2 + 9aP + 96P 2 , 



,42 

P 



where = p? + A is the effective mass of the scalar field and G^{t\,t2\p) is the full 
momentum space two point correlation function of the scalar field to leading order in 1/N. 
Fields evaluated at the saddle point are denoted by a bar. 

Note that G^(t±,t2;p) depends on the effective mass m^, and therefore it is difficult 
to solve (2.9) in full generality. Hence, in what follows we use the approximation proposed 
in [13]. In particular, we assume that tends to a stationary value m*^ and that this 
happens fast enough to be approximated by a jump. Then the two point correlation 
function G</,(ii,i2iP) is approximately the same as the propagator in the massive free field 
theory in which the physical mass is instantaneously changed from p,Q to m*±. i.e., 



G <j> (ti,t 2 ;p) ^ G^(ii,t 2 ;p;^o,m£) , (2.10) 



where [13] 



G^tuhw^ml) = (t f /° p)2 cos^(t 1 -t 2 )+ C f^ 2 ^ coso^x+^H-^e-^-fal , 

(2.11) 



with ujp = \Jp 2 + fn*^ and ojq p = \J p 2 + p^. The second term on the right hand side is the 
only one that breaks time translation invariance. However, its contribution to p vanishes 
for t\ = t 2 = t — > 00 within our stationary phase approximation. Therefore (2.9) yields the 
following equation for m*^ 

*2 2 94 f 1 r~T 2 2 ( W \ . \ 

m * = * ~ 4^ + 2 V * ~ K arccos [^J ~ A ) 

+ T& ("> + \yl m ? ~ & arccos (4) " A ) 2 ' (2 - 12) 

where we have taken a sharp cut off A to regulate the divergent integral over the momentum. 
To eliminate the cut off dependence we apply the following renormalization scheme 

2 2 1 94 a 96 .2 
= » + 4^ A+ 16^ A ' 

9? = 94 + £A . (2.13) 
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As a result, the gap equation for m*^ becomes 

m * 2 = Ml - J (iMi + \\j m f ~ Po araxwfrto/mj)^ 

(/Mi + \\j m f- Po arccos(^o/mj)^ . (2.14) 



+ 96 



16vr 2 

Solutions of this gap equation describe the stationary points of the effective potential. 
In the following we analyze these solutions and demonstrate that the quenched model 
exhibits peculiar phase structures. 

3 Phase structure of the model 

To analyze the admissible phases of the model let us derive the effective potential of the 
theory at t — > oo. From (2.8), we get, up to A-independent constant, 

v effiP, m <t> ) = Y p T P ~6 P 2 2 J J ?27r)2 G ^(*'*;^/ / 0' m ) • t 3 - 1 ) 

Varying this effective potential with respect to p and m^ 2 correctly reproduces the saddle 
point equations (2.9). Note that p is not a dynamical field since it enters only algebraically 
into the action (2.8). Hence we eliminate it from the effective potential using the second 
equation (2.9). Replacing the couplings by renormalized ones and further rescaling them 
by no yields 



V e ff(m 2 ) = — ^^^1 + — \/ m 2 — 1 arccos(l/V m 2 )^ 



~9f 


4(4vr; 


) 2 


96 




6(4tt; 






+ 



+ jtx3^ ( 1 + - \/~m 2 — 1 arccos(l/v / m 2 ) 



2 



3 



1 + — v m 2 — 1 arccos(l/v / m 2 ) J 
(m 2 + 2)y/m 2 — 1 arccos(l/V?Ti 2 ) + m 2 — logm 2 



48tt 



, (3.2) 



where V e ff, pn, gf and m 2 denote respectively the rescaled dimensionless effective poten- 
tial, the dimensionless renormalized couplings and asymptotic mass. 

Let us briefly discuss the case where both g§ and pr are zero. This case was considered 
in [13]. The characteristic shape of the effective potential of the quenched theory in this 
case is shown in figure l(red line), and as pointed out in [13] a finite mass always emerges 
at late times in the presence of interactions. This should be contrasted with the presence 
of a global minimum at m = in the unquenched theory as shown in figure l(blue line). 
This case already illustrates the main point that we wish to make, namely that the shape 
of the effective potential at late times depends on the quench, an event that occurred in the 
far past. More interesting and spectacular however is the case when the theory is sitting at 
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Figure 1: Red line (bottom): effective potential (3.2) as a function of m for g§ = [ir = 
and gf = 1. 

Blue line (top): effective potential of the t/> 4 theory, i.e., ge = 0, in the absence of quench 
when fiR = and g± = 1. 



the tricritical point, i.e., when fj,R = g^ = 0. Expanding effective potential then for large 
and small values of m? yields 



where g c = 256 corresponds to a critical value beyond which the potential is unbounded 
from below and thus the theory is unstable. 

It is remarkable that g c is larger than the corresponding value in the unquenched case 
[18]. There, the region of stability is bounded 3 by < g% < (4-7r) 2 = g§. In particular, 
previous studies[15, 16, 18] have spelled disaster for the theory in the ultraviolet limit: the 
/3-function drives the system into a UV fixed point g§ = 192 which lies beyond the region 
of stability. In contrast, our results indicate that there is a way to circumvent the above 
conclusion by a quench in the parameters of the system. 

It is instructive to contrast the phase diagrams with those of the unquenched case. 
If the coupling constant belongs to the range < g§ < (4-7r) 2 and the system is not 
quenched, then there is only one admissible conformal phase. Quenching the system in 
this regime, we explicitly break the conformal invariance by introducing a scale hq, and as 
a result the system resides in the light phase - a unique vacuum state associated with the 
minimum of effective potential, see solid blue graph at the top of figure 2 which represents 
a characteristic plot of the effective potential in this case. 

On the other hand, if the coupling constant is tuned to the special value g^ the un- 
quenched potential becomes flat, and thus a continuum of massive solutions emerges. This 

3 The lower bound g% > is necessary to avoid classical instability. 




if m >> 1 



if m « 1 . 



(3.3) 



-6- 




continuum is associated with spontaneous breaking of scale invariance and it coexists with 
the conformal (or massless) phase which we analyzed before. Quenching the system breaks 
the scale invariance explicitly and fio singles out a unique vacuum out of the continuum 
of states. As a result, in the steady state as t — > oo we find two admissible vacua shown 
in figure 2. The heavier phase for this value of the coupling is only meta-stable unless 
the coupling constant is sufficiently large, see dashed green and solid black graphs at the 
bottom of figure 2. 

Of course, in the unquenched case all phases become unstable for g§ > g§ and the 
system rolls down to infinity in a finite time [20]. Here we have showed how the quantum 
quench may enhance the stability of the system such that for g§ > g§ these phases are 
stable and the escape to infinity is avoided until g§ hits g c . 

4 Concluding remarks. 

To conclude, by studying the late time phase structure of the </> 6 theory after a quantum 
quench, we have essentially demonstrated in a specific example the following: a dramatic 
event that occurred in the far past can have significant effects even in the far future. Not 
entirely unexpectedly, we find that in the large N limit the late time physics cannot be 
described by simple thermalization with a single effective temperature, which has been 
noted in many integrable models. Explicit details are relegated to Appendix A. Instead of 
simply thermalizing, the quantum quench appears to modify the phase structure, even long 
after the system has relaxed and settled into an equilibrium state. From the calculations, it 
appears that this is a generic feature of quantum quenches, and is not specific to the model 
that we have studied. In passing, we note also significant modifications in a supersymmetric 
version of this model: contrary to the current model, where new stable phases are created, 
we found instability generated by the quench [21]. This dependence of the phase structure 
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on past events may have implications in other areas of physics, e.g. in cosmology. It is 
therefore important to determine and perhaps classify, different time dependent changes in 
a generic theory that could potentially lead to drastic modification of late time physics.. In 
this work we have made extensive use of techniques developed in [13], where it is implicitly 
assumed that the system relaxes ultimately to an equilibrium. The authors of [13] support 
their claim by explicit numerical computations and show that their assumption works very 
well in cj) theories in arbitrary dimensions. In this work we extrapolate their assumption 
to study the c/> 6 theory at the tricritical point. Further numerical checks of the current 
model and more general ones are under way and will appear elsewhere [21]. 

A Absence of thermalization: 

A question regarding the late time behavior studied in this work is whether it is effectively 
thermal. However, despite carrying interaction terms, the model we are studying is in fact 
integrable in the large N approximation, and as noted in previous works, is not expected 
to be describable by a single effective temperature [22]. In the rest of this appendix we 
demonstrate that the late time physics of our large N vector model is incompatible with 
simple thermalization leading to a single effective temperature. 

If the stationary behavior of the system is thermal, i.e., the effective mass is the 
thermal mass of the system. One can fix the temperature by matching the gap equation 
(2.9) at t — > oo with that in the thermal case. As a result, one gets the following relation 
which determines the inverse temperature (3 

(<p 2 )o = (^ 2 )/3 , (A.l) 

where subscripts "/3" and "0" indicate that the expectation values are taken in the thermal 
and l^o) states respectively. 

Having fixed the temperature via (A.l) (implicitly), we would like to compare the 
expectation value of the stress tensor of the thermal state with our quenched state at 
late times. The energy-momentum tensor of the scalar O(N) vector model in Minkowski 
signature is given by 

V = drf ■ d v cj> - r/^Qc^ • d»<j> - NV(tf/N)^j - - 7^,3 V . (A.2) 

where V((j) 2 /N) is the potential 4 . The Euclidean form T^ u is obtained by flipping the sign 
of the potential and replacing rjau with 5^ u in the above expression. In the large N limit 
the expectation value (V((p 2 /N)} equals V {(4> 2 ) / N) . Hence, provided that (A.l) holds, we 
get in the limit t — > oo 

(Too)o + (T e )t = (rioaVop - \i}ap\ W • d^) + (<W% - \s a ^j (d a <t> ■ d^) p . (A.3) 

The left hand side of the above expression should vanish if the system thermalizes since 
the zero-zero component of the Euclidean energy-momentum tensor is the minus energy 



4 In our case V^/N) = ^ + ^ 2 ) 2 + ^(0 2 ) 3 
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Figure 3: The plot of the free energy density / for various values of the coupling constant 
g§ as a function of effective mass m*i. N = 1, the inverse temperature is set to (3 = 1 and 
g% = 0.9(?g, #g, 1.05gg, l.lgg from top to bottom graph respectively. 

density. Consider adding rn^cj) 2 /2 to the first term on the right hand side. We immediately 
identify that as (Too)o of a quenched free scalar field whose mass jumps from to m\. On 
the other hand, subtracting m^(j) 2 /2 from the second term gives the minus thermal energy 
of the free scalar field with mass m*i. Since the free theory does not thermalize, there is no 
temperature such that these two terms cancel each other. This proves a mismatch between 
thermal and quenched physics. 

It is also instructive to compare the thermal free energy of the system at the tricritical 
point with the corresponding effective potential (3.2). In the large N limit the thermal free 
energy density / can be evaluated in a closed form and is given by 

f/ N =9a( 1 ^Z + log(l-e-^) \ 3 _ n4 _ m 2 T \og{\ - e-^) 

Jo 47r P 

The typical plots are shown in figure 3. While for g$ < g§ the free energy exhibits cer- 
tain similarity with the corresponding effective potential, they become manifestly different 
for g$ > and hence thermalization is not expected. 
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